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We provide a detailed description of our previously proposed scheme for topological interface engi- 
neering with constructed defects and textures perforating across coherent interfaces between different 
broken symmetries [M. O. Borgh and J. Ruostekoski, Phys. Rev. Lett. 109, 015302 (2012)]. We 
consider a spin-1 Bose-Einstein condensate, in which polar and ferromagnetic phases are prepared 
in spatially separated regions. We show that a stable coherent interface is established between the 
two phases, allowing defects of different topology to connect continuously across the boundary. We 
provide analytic constructions of interface- crossing defect solutions that could be experimentally 
phase-imprinted using existing technology. By numerically minimizing the energy, we calculate the 
core structures of interface-crossing defect configurations. We demonstrate nontrivial core deforma- 
tions to considerably more complex structures, such as the formation of an arch-shaped half-quantum 
line defect, an Alice arch, at the interface, with the topological charge of a point defect, whose emer- 
gence may be understood by the "hairy ball" theorem. Another example of an energetically stable 
object is the connection of a coreless vortex to a pair of half-quantum vortices. We show that rota- 
tion leads to spontaneous nucleation of defects in which a coreless vortex continuously transforms 
to a half-quantum vortex across the interface. 
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I. INTRODUCTION 

In systems described by an order parameter, for ex- 
ample superfluid liquid helium, liquid crystals and Bose- 
Einstein condensates (BECs), the symmetry properties 
of this order parameter determine the topological proper- 
ties of defects and textures [1]. Systems with vector order 
parameters, such as superfluid liquid ^He [2] and spinor 
BECs [3|-l5|, may exist in several distinct phases. Each 
phase corresponds to a different ground-state manifold 
of energetically degenerate and physically distinguishable 
states, resulting in different order-parameter symmetries. 
When two phases coexist in spatially separated regions in 
the same medium, topological defects cannot penetrate 
the boundary unchanged, but must either terminate at 
the interface or connect nontrivially to an object on the 
other side. Topological interfaces appear as important 
objects in seemingly distant areas of physics. For ex- 
ample, it has been proposed that a series of symmetry 
breakings in the early universe leads to the formation of 
cosmic strings that terminate on boundaries between re- 
gions of different vacua [gI, . Complex interface physics 
also arises in string theory from collisions between branes 
during inflation [S^, ^] , in condensed-matter theory in ex- 
otic superconductivity [lo|, and in superfluid liquid ^He 
when a magnetic-field gradient causes A and B phases to 
coexist, resulting in the pos sibility of nontrivial defects at 
the phase boundary |lll-[l^ . Parallels between defects in 
superfluids and objects in cosmology [7] led to the sug- 
gestion of using superfluid systems to study analogues 



* M.0.Borgh@soton.ac.uk 
"f" j anne@soton. ac. uk^ 



of cosmo logical phenomena in the laboratory |12l, 15L 
including defect formation in phase transitions |16l-[l8|, 
analogues of brane annihilation [l^, ll9|] , and structures 
similar to cosmic vortons 

Here we present a detailed description of our proposal 
for an experimentally feasible scheme to study topolog- 
ical interface physics in a gas of ultracold atoms with 
spin degree of freedom [23]. A stable, coherent interface 
between two ground-state manifolds of different broken 
symmetries may be prepared by local manipulation of 
the interaction strengths of different scattering channels 
of the atoms. The proposed technique can be used to 
realize, e.g., two different ground-state manifolds simul- 
taneously in spatially separate regions. Defects and tex- 
tures crossing the interface, and connecting objects of 
different topologies, can be created by controlled phase- 
imprinting of elementary vortex and soliton structures 
or by rotating the gas. Under dissipation the phase- 
imprinted defect configurations relax to more complex 
objects for which the continuous spinor wave function 
interpolates smoothly across the topological interface. 

As an example we consider defects and textures cross- 
ing the boundary between polar and ferromagnetic (FM) 
regions in a spin-1 BEC and show that a coherent in- 
terface is established within a continuous condensate 
wave function. We analytically construct prototype 
spinor wave functions representing interface-crossing de- 
fect structures, and by numerically minimizing their en- 
ergy, we evaluate the configurations that emerge as a 
result of energy dissipation in spin-1 BECs. The sim- 
ulations demonstrate nontrivial core deformations and 
defect structures. We characterize the defect cores, ana- 
lyze the energetic stability of defect solutions crossing the 
interface, and explain how the defect-carrying conden- 
sate wave function continuously interpolates across the 
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interface. We also demonstrate nucleation of interface- 
crossing defects consisting of a coreless vortex that con- 
nects to a half-quantum vortex. 

A spinor BEC |3M5| is created in an all-optical trap so 
that the spin degree of freedom of the atoms is not frozen 
out by magnetic fields. Spin rotations then combine with 
the condensate phase to form a large set of physically 
distinguishable states. Because also the contact interac- 
tion between the atoms in the condensate becomes spin- 
dependent, energetically degenerate subsets depend on 
the strength and sign of the spin-dependent contribu- 
tions. Hence the spinor BEC exhibits a rich diagram 
of phases with different broken order-parameter symme- 
tries (23-[3Qj as a function of the interaction strengths. 
This is similar to superfluid liquid ^He where nonzero 
spin and orbital angular momenta of the Cooper pairs 
combine to form phases with different order-parameter 
symmetries supporting a variety of defects and tex- 
tures [31^. 

Modern techniques used in experiments with ultracold 
atoms provide tools for unprecedented control over sys- 
tem parameters and for accurate measurements, includ- 
ing the possibility for in situ observation of vortices in 
spinor BECs. As there has been considerable interest in 
the studies of the stability properties of field-theoretical 
solitons in various physical systems (32l-[35| , it is therefore 
not surprising that this is also followed by an accelerating 
theoretical interest in a variety of stable and meta-stable 
objects in multicomponent atomic BECs. Perhaps the 
simplest of such structures where the multicomponent 
nature of BECs plays an important role are ID vector 
solitons (36l-[43|, typically consisting of stable combina- 
tions of dark and bright solitons in different condensate 
components. Higher dimensional defects and textures in- 
clude vortex sheets [44] and 3D particle-like solitons [20- 
[22, 45-47] in two-component (pseudo-spin- 1/2) conden- 
sates, as well as a rich phenomenology of defects and 
textures in spin-l m, m, m, , spin-2 [H-ii] , and 

spin-3 [30|, E3 BECs. Interface physics has been stud- 
ied in the two-component BEC systems, for example in 
the context of vortex bifurcation at energetically estab- 
lished interfaces in the phase-separation regime [68|, [69| 
and interface collisions [19]. There is a rapid parallel 
experimental development, exemplified by preparation 
of coreless vortices and related textures |7Q|-[72]. as well 
as observations of singular vortices produced in phase 
transitions [73], and of spin-texture formation [74-76]. 
Furthermore, trapping of ultracold atoms in artificial 
gauge-field potentials [73 was recently realized experi- 
mentally [78, 79]. This presents intriguing possibilities 
for the stability studies of defects and textures, including 
those of particle-like solitons [47^. 

Our study of topological interface engineering is orga- 
nized as follows: In Sec. [Ill we first give a brief overview 
of the standard mean- field theory of the spin-l BEC and 
then proceed to give a more detailed presentation of the 
topology and basic defects of the FM and polar phases. 
In Sec. mil we discuss how an interface between FM and 



polar regions can be created, and then identify and ex- 
plicitly construct interface-crossing defect solutions. We 
proceed to minimize the energy of the defect solutions in 
Sec.HVland describe the emerging structures of the defect 
core and the energetic stability of the defects. We explic- 
itly demonstrate continuity of the spinor wave function 
across the stable interface. We summarize our findings 
in Sec. El 



II. SPIN-1 BEC 

A. Mean-field theory of spin-l BEC 

We consider a trapped spin-l atomic BEC confined 
in an all-optical, harmonic trap. We may then employ 
the classic Gross-Pitaevskii mean-field theory describing 
a spatially inhomogeneous macroscopic condensate wave 
function ^(r). Since we are considering spin-l atoms, 
^(r) can be written in terms of the density of atoms 
n(r) and a normalized, three component spinor C(r) in 
the basis of spin projection onto the z axis as 



^(r) = Vn(r)C(r) = vM^) 




Ck = i. (1) 



The mean-field Hamiltonian density then reads [23, [25|, 
80] 
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where V{t) is an external trapping potential and m is the 
atomic mass. In this work we consider the atoms trapped 
in a slightly elongated potential, so that 



V(r) = -muj^ I X 
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The spin operator F is given by a vector of spin-l Pauli 
matrices. Its expectation value (F) = Ci^a^C^ is the 
local spin vector. A weak external magnetic field may 
be imposed, in which case linear and quadratic Zeeman 
shifts as described by the last two terms will arise. Most 
of our numerical results correspond to the cases for which 
the Zeeman splitting is assumed to be negligible. We 
note, however, that all our results remain qualitatively 
the same in the presence of weak Zeeman splitting energy. 

We also investigate the configurations of defects and 
textures in a rotating trap. In that case we minimize the 
free energy in a rotating frame, corresponding to 



H' = H - {rt'V) 
= Id'r 
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where we have assumed the axis of rotation defined by ft 
to be along the positive z axis, and L denotes the angular 
momentum. 

The two interaction terms in Eq. (|2]) arise from the fact 
that the spins of two colliding spin-1 atoms may combine 
to either or 2. There are therefore two s-wave scattering 
channels, with scattering lengths ao and a2, contributing 
to the contact interaction between the atoms in the con- 
densate. Standard angular- momentum algebra Q sep- 
arates the interaction energy into one spin-independent 
contribution and one term that depends on the magni- 
tude of the spin. The strengths of the spin-independent 
and spin-dependent interaction terms are then given by 
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respectively. Additional magnetic dipole-dipole interac- 
tions that may influence the spin textures jH, [ll, [HI are 
neglected here. 

The sign of C2, the strength of the spin-dependent 
interaction, determines the magnitude of the spin vec- 
tor in a uniform ground state, leading to the two topo- 
logically distinct phases of the spin-1 BEC. If C2 < 0, 
energy minimization favors maximized spin magnitude 
|(F)| = 1 in the FM phase. This is the case for ^'^Rb 
where co/c2 — —216 [82|- Conversely, if C2 > 0, as for 
2^Na with co/c2 :^ 31 [sl, |(F)| = is favored in the 
polar phase. The two phases are described by funda- 
mentally different order parameters, supporting different 
families of defects, which we will discuss in some detail 
below. 

Characteristic length scales arise from the interac- 
tion terms. The spin-independent interaction defines the 
usual density healing length 
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which describes the length scale over which the atom den- 
sity n(r) heals around a local density depletion. In addi- 
tion, the spin-dependent interaction gives rise to a spin 
healing length 
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defining the distance over which |(F(r))| heals as the or- 
der parameter is excited out of its ground-state manifold. 
That situation arises in two cases of importance for the 
analysis presented in this article. Firstly, the core of a 
singular vortex in one phase may fill with atoms such that 
the atoms at the singularity exhibit the opposite phase. 
This can happen since the singularity of the spinor order 
parameter may be accommodated either by forcing the 
density to zero or by requiring that the wave function 
become orthogonal to the ground-state manifold (mean- 
ing locally perturbing |(F)|) [56, 63]. The size of the 
filled vortex core is then determined by ^f- Secondly, 



we are interested here in the interface between polar and 
FM regions. For the condensate wave function to inter- 
polate between the two manifolds, the spin magnitude 
must leave its ground-state value close to the interface, 
which will therefore acquire a width determined by ^f- 



B. Ground state manifolds and basic defects 

The order parameter manifold is the set of energeti- 
cally degenerate, physically distinguishable states. In the 
condensation transition, this symmetry is spontaneously 
broken, and this broken symmetry determines the topo- 
logically distinct families of defects. The FM and polar 
phases of the spin-1 BEC are described by very differ- 
ent order-parameter manifolds, leading to dramatically 
different possible vortex states. Before discussing the in- 
terface between FM and polar regions in the next section, 
we here give an overview of the families of defects in the 
purely FM or purely polar BEC. 



1. The FM phase 

If C2 < in Eq. (j2j), the spin-dependent interaction 
will favor a state that maximizes the magnitude of the 
spin everywhere, such that |(F)| = 1. A representative 
FM spinor is given hy ( = (1,0,0)^, such that the spin 
vector is parallel with the z axis. From this representative 
spinor, a general FM spinor may be constructed by a 3D 
spin rotation 

U (a, /3, 7) = exp(-iF^a) exp(-zF^/3) exp(-zF^7) , (8) 

defined by three Euler angles, together with a conden- 
sates phase ^, as 



*^/7(a,/3,7) 



' cos^ f 
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V2e'^ sin^ I 



(9) 

where the condensate phase is absorbed in the third 
Euler angle: 7' = 7 — (/). Any FM spinor is thus de- 
scribed by some particular choice for (a,/3,7'). Thus 
the broken symmetry of the ground-state manifold 
is represented by the group of 3D rotations SO (3). 
The spin vector is given by the Euler angles as 
(F) = (cos a sin /3, sin a sin /3, cos /3) . 

Topological stability of line defects can be character- 
ized by studying closed contours around the defect line 
and the mapping of these contours into order-parameter 
space [11]. If the image in order-parameter space of a 
closed loop encircling a line defect can be contracted to 
a point, the defect is not topologically stable. The order- 
parameter space of a FM spin-1 BEC, SO (3), can be 
represented geometrically as (the unit sphere in four 
dimensions) with diametrically opposite points identified. 
If a closed contour connects such identified points more 
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than once, any pair of such connections can be ehminated 
by continuous deformation of the contour. Therefore any 
contour with an even number of connections can be con- 
tracted to a point, whereas a contour with an odd number 
of connections can be deformed into a contour with just 
one connection. Hence we have only two distinct classes 
of vortices: singly quantized, singular vortices that corre- 
spond to noncontractible loops, and nonsingular, coreless 
vortices representing contractible loops. All other vor- 
tices can be transformed to either one of these by local 
deformations of the order parameter. Mathematically, 
these equivalence classes are characterized by the two el- 
ements of the first homotopy group, 7ri[SO(3)] = Z2. 

The simplest representative of the class of singular line 
defects is constructed as a 27r winding of the condensate 
phase such that = —cp in Eq. ([9]), where cp is the az- 
imuthal angle. This results in a line singularity oriented 
along the z axis. The spin texture is uniform, such that 
a and P are arbitrary but constant. The vortex is then 
described by the spinor 
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From other vortices in the same equivalence class 
can be formed by local spin rotations. For example, we 
may rotate the spin vector such that at each point it 
points radially away from the vortex line, as illustrated in 
Fig. m^a). This vortex corresponds to the choices a = (p^ 
P = 7r/2 and 7' = in Eq. (|9]), and is described by 
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This singular spin vortex [2J, |25| illustrates another im- 
portant aspect of the FM phase: circulation alone is not 
quantized. The superfiuid velocity in the FM phase 



V = ( V7' + cos /3Va) , 

m 



(12) 



vanishes when 7' = and /3 = 7r/2. Therefore the circu- 
lation is in fact zero in C®^, but it does carry a nonvan- 
ishing spin current around a singularity of the FM spin 
vector whose structure is similar to an analogous vor- 
tex with a radial disgyration of the angular momentum 
vector in ^He. Further local spin rotations yield other 
singular vortices with different spin structures, such as 
the cross disgyration shown in Fig. [Hb) or a tangential 
disgyration with (F) = cp. 

A striking manifestation of the nonquantization of cir- 
culation in the FM phase is the formation of a non- 
singular coreless vortex. This can be constructed as a 
combined rotation of the spin vector and the condensate 
phase [Fig. Hl^c)]: a = —7' = (p, yielding the spinor 



V2cos^ ^ 



(13) 



where the Euler angle /3 is now a function of the radial 
distance p = \/ + such that /3 ^ as p ^ 0, keep- 
ing the spin texture continuous. The superfiuid velocity, 
Eq. (p!2]) . becomes 



mp 



cos / 
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and increases smoothly from zero at p = as /3 increases 
away from the vortex, the spin vector forming a fountain- 
like texture. The coreless vortex in the FM phase of a 
spin-1 BEC is analogous to the Anderson- Toulouse and 
Mermin-Ho vortices in superfiuid ^He [13, [HI, which dif- 
fer by the boundary condition imposed on the angular 
momentum vector at the container wall. In the BEC 
there is no hard container wall, and the amount by which 
P turns from the vortex line to the edge of the cloud is 
determined by the rotation of the trap, causing the total 
angular momentum to vary smoothly with rotation. 

The coreless vortex can be continuously transformed 
into other members of the class of nonsingular vortices, 
including the vortex- free state, by purely local opera- 
tions. The continuous deformation is a striking conse- 
quence of the two-element character of the fundamental 
group of the SO (3) order-parameter space: the doubly 
quantized vortex belongs to the same topological class as 
the nonsingular vortices and the vortex- free state, and 
can be continuously unwound, if the orientation of the 
spin texture is not fixed outside the structure. Another 
nontrivial nonsingular vortex with continuous spin tex- 
tures is displayed in Fig. [TJd). 

We have now identified two topologically distinct 
classes of vortices that can both carry mass and spin 
circulation in the FM spin-1 BEC. A similar situation 
applies in the A phase of ^He. There the singular vortex 
has the lower energy, but the energy barrier for nucleation 
is lower for the nonsingular vortex [86| . In a rotating FM 
spin-1 BEC, the coreless vortex has the lower energy and 
the lower nucleation barrier, and consequently the ground 
state is made up of non sing ular coreless vortices for suf- 
ficiently rapid rotation |53i 153, IHsl , IHol , [6l|. However, it 
is also possible to form a singly quantized, singular vor- 
tex [51, 52], which despite not being the lowest-energy 
state at any frequency of rotation can nevertheless be 
energetically stable as a local energy minimum [63]. 

So far we have considered line defects, classified by the 
first homotopy group tti. Point defects — monopoles — 
are analogously classified by the second homotopy group 
7r2. For the FM order-parameter manifold, the second 
homotopy group is the trivial group, 7r2[SO(3)] = 0, in- 
dicating that the FM phase does not strictly speaking 
support point defects. However, it is possible to form a 
spinor with a monopole structure of the spin vector (a 
radial hedgehog) as the termination of a doubly quan- 
tized vortex [53 (Fig- [2|). This is the analogue of the 
Dirac monopole in quantum field theory, and the dou- 
bly quantized vortex line is called the Dirac string. The 
corresponding spinor is written by choosing a = ^' = (p 
and /3 = 6 (where 6 and (p are the polar and azimuthal 
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FIG. 1. (Color online) Nontrivial vortices in the FM phase, 
(a) A radial disgyration of the spin vector (red cones) around 
the singular core represents a spin vortex. This vortex carries 
a spin current but no mass circulation, a manifestation of the 
nonquantization of circulation alone in the FM phase, (b) 
All vortices in the class of singly quantized, singular vortices 
can be deformed into each other by local spin rotations. For 
example, it is possible to form a spin vortex with a cross- 
disgyration of the spin vector, (c) In addition to the one class 
of singular vortices, the FM phase also supports nontrivial, 
nonsingular vortices. In a coreless vortex the order parame- 
ter remains nonsingular everywhere in the vortex structure. 
The characteristic fountain structure of the spin is formed by 
a rotation of the spin vector around the vortex line, together 
with a winding of the condensate phase, corresponding to spin 
rotations about the local spin vector (indicated by the orthog- 
onal green and blue vectors), (d) Like the class of singular 
line defects, members of the family of nonsingular vortices are 
related by local deformations of the order parameter allowing 
different vortex configurations. 



angles, respectively) to form 
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This Dirac monopole can be continuously deformed into 
the spin structure of the coreless vortex [57,] . 



2. Polar phase 

We next consider C2 > in Eq. ([2]), in which case the 
spin-dependent interaction favors a state with |(F)| = 
0. A simple representative polar phase spinor fulfilling 
this requirement is = (0, 1,0)^. As for the FM phase, 
the general polar spinor is found by applying a 3D spin 
rotation U (a, /3, 7) together with a condensate phase (j): 



V2 



-e sin P 
V2COS/3 
e*" sin P 



(16) 



FIG. 2. (Color online) An analogue of the field-theoretical 
Dirac monopole can be constructed in the FM phase of the 
spin-1 BEG. The radial hedgehog monopole structure of the 
spin vector (cones) appears as the termination of a doubly 
quantized, singular vortex line. 



We now make the important observation that the unit 
vector d = (cos a sin /3, sin a sin /3, cos /3) defines the local 
direction of macroscopic condensate spin quantization. 
This allows us to rewrite in terms of this vector as j56[ 



V2 



d'jQ I xdy 

dx ^dy 



(17) 



The condensate phase ^, which takes values on a unit 
circle, and the unit vector d, taking values on a sphere, 
thus together fully specify the order parameter in the po- 
lar phase. Note, however, that C^(^, d) = C^(0 + 7r, — d). 
These two states must be identified in order to avoid 
double counting. The order parameter space therefore 
becomes (U(l) x S'^)/Z2, where the factorization by the 
two-element group Z2 results from the identification. The 
vector d should thus be taken to be unoriented and de- 
fines a nematic axis [2^] , and the order-parameter is cor- 
respondingly said to exhibit nematic order, which leads 
to parallels with the A-phase of superfluid ^He. 

A simple singly quantized vortex can again be con- 
structed as a 27r winding of the condensate phase, keeping 
d uniform (choosing a and P to be constants): 



-e sin p 
= ^ I V2cosp 
e*" sin P 



V2 



(18) 



(19) 



In the polar phase the superfluid velocity is 

V — V^. 

m 

We observe that v depends only on the gradient of the 
condensate phase, and is independent of d. This means 
that another singly quantized vortex, with the same cir- 
culation as that described by Eq. (p!8|) . can be formed 
by allowing d to wind by 27r (thus preserving single- 
valuedness of the order parameter) in addition to the 
winding of the condensate phase. This is achieved by 
choosing a = (j) = (f in Eq. (p!6|) . yielding the spinor 



— sin P 
V2e^^ cos/3 
e^*^ sin/3 



(20) 
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FIG. 3. (Color online) The polar phase of the spin-1 BEC 
supports half-quantum vortices, constructed as a tt winding 
of the condensate phase (which determines the quantized cir- 
culation) together with a d ^ — d winding of the nematic 
axis. The identification ("P((/), d) = ^^(0 + tt, — d) keeps the 
order parameter single- valued, (a) Half-quantum vortex as 
described by Eq. ()2ip . The nematic axis (red cones) winds 
by TT in the plane perpendicular to the vortex Une as the core 
is encircled. The disclination plane where the d ^ — d iden- 
tification is made is indicated by the dashed line, (b) The 
winding of the nematic axis need not stay in the plane per- 
pendicular to the vortex line. For example, a half-quantum 
vortex can also be constructed where d winds in a plane par- 
allel to the vortex line as the defect is encircled. 

One can further show from Eq ([19]) that circulation is 
quantized in the polar phase. However, due to the ne- 
matic order, the smallest circulation possible is half that 
of a singly quantized vortex. The equivalence (^, d) = 
(■^(^ + TT, — d) implies that we can allow the condensate 
phase to wind by tt along a loop encircling the vortex and 
still preserve single- valuedness of the the order parame- 
ter by a simultaneous d ^ — d winding of the nematic 
axis [49]. If d is in the (x, ?/)-plane, a half-quantum vor- 
tex can be written 

^hq ^ e ^ = ^ . (21) 

In general, the axis about which d winds need not co- 
incide with the vortex core. Figure ^h) shows a half- 
quantum vortex where d winds about an axis perpen- 
dicular to the vortex line. This vortex is related to that 
shown in Fig [3fa) and defined by Eq. (|2T]) by a spin ro- 
tation. The resulting spinor wave function may appear 
quite complicated, but the tt winding of the nematic axis 
still allows us to identify the vortex. 

Thus circulation is quantized in the polar phase, and 
indeed one can show that 7ri[(U(l) x S'^)/Z2] = Z. The 
topological charges 1/2, 1, 3/2 etc. are additive. For 
example, the state with two half-quantum vortices be- 
longs to the same equivalence class as the singly quan- 
tized vortices. This observation shall prove important for 
understanding the core structure of defects that cross a 
polar-FM interface. 

In addition to singular line defects, the polar phase 
also supports singular point defects [50|, Ell: A spheri- 
cally symmetric point defect, analogous to the 't Hooft- 
Polyakov monopole in quantum field theory, is formed 



by choosing d = f = (sin ^ cos sin ^ sin cos 6>) in 
Eq. (p!6|) , such that the vector field d forms a radial hedge- 
hog structure 

/ -e-^'^ sin l9 \ 
^pm= V^cos^ . (22) 

y e^^sinO J 

Here a singular point defect is located at the origin. The 
two spinor wave function components (± form overlap- 
ping, singly quantized vortex lines with opposite circu- 
lation. The vortex lines are oriented along the z axis, 
normal to a dark soliton plane (phase kink) in the com- 
ponent Co- 



lli. TOPOLOGICAL INTERFACE IN A SPIN-1 
BEC 

The two phases of the spin-1 BEC exhibit different 
topological properties of the order parameter, which is 
manifest in the very different defects they support, as 
detailed in the previous section. We now consider the 
behavior of the order parameter and defects when the two 
phases are realized simultaneously in spatially separated 
regions within the same ultracold gas, so that an interface 
between the different ground-state manifolds is formed. 

The order-parameter manifolds in the FM and polar 
phases emerge out of the full symmetry of the condensate 
wave function as the spin-dependent interaction selects 
an energetically degenerate subset of all possible wave 
functions. In particular, the order-parameter space in the 
FM phase consists of all spinor wave functions that max- 
imize the spin magnitude, |(F)| = 1, everywhere if the 
texture is uniform. Correspondingly in the polar phase, 
the order-parameter space is the set of wave functions 
that have |(F)| = everywhere. These sets are clearly 
nonover lapping. However, because they form subsets of 
the same spin-1 wave function, a continuous connection 
between spatially separated polar and FM manifolds is 
possible by exciting the wave function out of its ground- 
state manifold close to the interface, locally restoring its 
full symmetry. 

In this section we first discuss how the interface may 
be created in the spinor BEC through local manipula- 
tion of the scattering lengths. We then identify the basic 
interface-crossing defect states and explicitly construct 
prototype spinor wave functions to describe them. In 
Sec. lIVI we minimize the energy of these spinor wave func- 
tions for defect configurations and show how this leads 
to a rich phenomenology of defect structures. 



A. Creating the topological interface 

In order to realize a topological interface in a spinor 
BEC it is necessary make the system switch continuously 
between regions of different broken symmetries. Which 
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broken symnietry is preferred in a spin-1 BEC on ener- 
getic grounds is determined, as explained in Sec. Ill A| 
by the spin-dependent interaction. Therefore if one can 
spatially control the interaction strength C2, Eq. (|5j), sep- 
arate FM and polar regions within the same BEC can be 
engineered. Specifically, since C2 oc (a2 — ao), this implies 
changing the ratio ao / (22 of the two scattering lengths of 
colliding spin-1 atoms such that C2 changes sign. 

The scattering lengths that determine interaction 
strengths in ultracold-atom systems are routinely manip- 
ulated using magnetic Feshbach resonances. However, 
this technique cannot be used for our present purpose 
since the strong magnetic fields required would freeze out 
the spin degree of freedom and destroy the spinor nature 
of the BEC. The possibility for engineering the scatter- 
ing lengths in the spinor BEC is instead provided by the 
use of eithe r op tical js^] or microwave-induced Feshbach 
resonances \8§\, in which case the fields can be kept suf- 
ficiently weak in order not to destroy the spinor nature 
of the BEC. 

We suggest constructing an interface between topolog- 
ically distinct manifolds in a spinor BEC by local adjust- 
ment of the scattering lengths, such that regions with 
different-sign C2 are created [23|. In a spin-1 BEC exper- 
iment, the spatial dependence of the scattering lengths 
can then result in an interface between coexisting FM and 
polar phases. Doing so, however, presents practical chal- 
lenges. If an optical Feshbach resonance is used to manip- 
ulate one or both scattering lengths, the spatial pattern 
corresponding to a sharp interface can be imposed using 
a holographic mask. Optical Feshbach resonances suf- 
fer from inelastic losses [83, but these can be kept small 
for small adjustments of the scattering lengths. Since 
the spin-dependent interaction is proportional to the dif- 
ference between ao and a2, only a small relative shift is 
needed to create the interface if |c2| is small, which is true 
for both ^^Rb and ^^Na, commonly used in spinor-BEC 
experiments. 

Using a microwave-induced Feshbach resonance avoids 
the problem of large inelastic losses, but makes engineer- 
ing the spatial profile more difficult, since, except in spe- 
cific traps, for example surface microtraps, a microwave 
field cannot be focused in the same way as the laser. 
This problem could be overcome by using an optically 
induced level shift to tune the microwave transition off- 
resonant. The microwave field could then be applied uni- 
formly across the system, whereas spatial control of the 
laser field is used to apply the optical tuning only in the 
region where no shift of C2 is required. 



B. Construction of prototype interface spinors 

The use of optical or microwave-induced Feshbach res- 
onances can thus realize spatially separated polar and 
FM phases in the same spin-1 BEC, with the condensate 
wave function remaining continuous, allowing, in prin- 
ciple, defects to connect across the interface. In order 



to demonstrate the nontrivial nature of defect penetra- 
tion across the interface between topologically distinct 
manifolds, we consider a spin-1 BEC where C2 abruptly 
changes sign at z = 0. We choose C2 > for z > and 
C2 < for z < 0, such that the interface exists at z = 
with the polar phase above it and the FM phase below. 
In the following, we analytically construct spinor solu- 
tions that represent physical wave functions for defects 
and textures simultaneously in the two different mani- 
folds. 

The simplest vortex connection can be identified by 
considering a singly quantized vortex in both phases, as 
illustrated schematically in Fig. SJ^a). Note that a singly 
quantized vortex does not mean the same thing in the 
two phases: the topology that describes vortices is en- 
tirely different, one vortex being a product of the broken 
symmetry manifold SO (3), with the fundamental homo- 
topy group of two elements, and the other one resulting 
from the broken symmetry (U(l) x 6'^)/Z2, with the fun- 
damental homotopy group of integers that represent the 
number of half-quanta of circulation. It is therefore not 
obvious that these different topological objects can be 
continuously joined across the interface. 

In the following we show how to construct spinor 
wave-function solutions that simultaneously represent a 
singly quantized vortex line in both phases and perforate 
through the interface with a 27r winding of the conden- 
sate phase around the vortex line (which we take to be 
along the z axis). A similar procedure is then extended 
to other topological defects and textures. The joining 
of two singly quantized vortex lines can be achieved by 
changing the sign of either of the spinor components ("+ or 
C_ . By appropriate choice of parameters doing so causes 
the spinor wave function to adjust between the two man- 
ifolds by forcing |(F)| to switch from to 1, or else leads 
to a state which immediately relaxes to the desired con- 
figuration. Physically, such a sign-change in one of the 
two spinor components can be obtained by introducing 
a dark soliton plane (phase kink) in that component at 
z = 0. The TT-phase shift across the soliton is then asso- 
ciated with a vanishing density in that spinor component 
at the soliton core. The BEC wave function, however, re- 
mains continuous across the interface, since the remain- 
ing spinor components have nonvanishing atom densi- 
ties also at the position of the soliton plane. The BEC 
wave function thus connects the two manifolds. In this 
construction, the switch between polar and FM sides is 
abrupt. In section IIVI we will see that as energy is re- 
laxed, the interface acquires a finite width determined 
by the spin heahng length Eq. ([7]). 

Following this procedure and starting from the expres- 
sion for a singular vortex in the FM phase, Eq. (p!Q]) . we 
can write the spinor wave function connecting two singly 
quantized vortices across the polar-FM interface explic- 
itly as 

C^-^ = — sin/3 , (23) 

^ VTV^e-sin^f y 
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FIG. 4. (Color online) Schematic illustrations of possible vor- 
tex connections. The polar phase is above the interface and 
the FM phase below, (a) The simplest vortex connection to 
construct is to consider a singly quantized vortex in either 
phase. Even though these are different objects, representing 
different topologies, their wave functions can be joined contin- 
uously across the interface, (b) If the singly quantized vortex 
in the polar phase also contains a 2it winding of the nematic 
axis, its spinor wave function can be made to connect to a 
spinor representing a coreless vortex in the FM phase, (c) 
The FM Dirac monopole has a spinor structure similar to the 
coreless vortex (in fact, the monopole can be continuously un- 
wound into a coreless vortex 57]) and it is therefore possible 
for it to form the termination point of a singly quantized po- 
lar vortex that ends at the interface, (d) Similarly, a polar 
monopole on the interface can form the termination point of 
a singly quantized FM spin vortex. The spinor is constructed 
by noticing that both defects have common phase windings 
of the spinor components. 



where the negative sign is used on the polar side of the 
interface and the positive sign on the FM side. Note that 
only the choice /3 = 7r/2 yields |(F)| = corresponding 
to an exactly polar state above the interface. However, 
even for a different j3 the spinor wave function has the 
appropriate vortex structure and will quickly relax to 
the polar phase for z > with a singly quantized vor- 
tex. This highlights the general consideration that even 
though writing exact vortex connections analytically may 
be very complicated, we have a simple method for find- 
ing approximate spinor wave functions representing de- 
fect connections. 

In Eq. (|2Q|) we demonstrated that a singly quantized 
vortex in the polar phase can include a 27r winding of 
the nematic axis in addition to the 2it winding of the 
condensate phase. Comparing Eq. (|2Q|) with the solution 
for a coreless FM vortex in Eq. ([13]), we note that these 
have a similar structure in terms of the complex phases 
of the spinor components. We can therefore construct an 
approximate wave function representing the connection 



of a singly quantized vortex in the polar phase with a 
coreless vortex on the FM side of the interface by the 
insertion of a soliton plane in in Eq. ([20]) . This yields 
the interface spinor 

/ -sin/3 \ 
^i^ci^ 726^^ COS /3 , (24) 

^2 y ±e2^^sin/3 ) 

where the positive sign is used on the polar side of the 
interface and the negative sign on the FM side. Choosing 
/3 = 7r/4 or /3 = 37r/4 yields |(F)| = 1 on the FM side, 
and specifically the choice {3 = 37r/4 approximates the 
coreless vortex [Eq. ([IS])]. This solution relaxes to the 
characteristic fountain-like spin profile to yield the state 
illustrated in Fig. |4]^b). 

In addition to connecting to another vortex across the 
interface, a vortex could also terminate on a point defect 
at the interface. Such solutions can be constructed by 
joining the monopole spinor wave functions of Eqs. ([15]) 
and ([22]) to vortices with analogous phase windings in 
each of the spinor components. For example, consider 
the Dirac monopole in the FM phase [Eq. ([T5]) ]. By the 
same construction that resulted in the interface-crossing 
defect in Eq. ([24]) we can connect the singular vortex of 
Eq. ([20]) on the polar side to the monopole of Eq. ([15]) 
at the interface by inserting a soliton plane into • The 
resulting spinor wave function, 

/ V2e-^'^ cos^ f \ 
C^^^ = — e-^^sin(9 , (25) 

V^V Tv/2sin^f / 

represents the monopole on the FM side of the interface. 
Here the negative sign refers to the polar side and the 
positive sign to the FM side. On the polar side the spinor 
has a structure similar to Eq. ([20]) , thus approximating 
a singly quantized vortex in the polar phase. The re- 
sulting defect configuration is illustrated in Fig. |4]^c) and 
represents a vortex in the polar phase terminating to a 
monopole at the interface. This defect is closely related 
to the one shown in Fig.jH^b), as the Dirac monopole and 
the coreless vortex can be deformed into each other by 
purely local operations. Note also that in Eq. ([25]) the 
Dirac string is represented by the singular polar vortex 
along the positive z axis and there is no line defect on 
the FM side. By instead aligning the Dirac string with 
the negative z axis, the doubly quantized line defect ter- 
minates on the monopole from the FM side, while for 
positive the spinor still represents a singular vortex 
connecting to the monopole from the polar side. 

In a similar way a point defect (radial hedgehog) in the 
polar side [Eq. ([22]) ] can be placed on the interface as the 
termination point of a singular FM vortex. We consider 
a defect structure with overlapping, singly quantized vor- 
tex lines in (± , both oriented normal to the interface and 
of opposite circulation, together with 7r-phase kinks in 
and Co- This spinor wave function can be parametrized 



as 

^sv^pm ^ V^COS^ , (26) 

using the negative sign on the polar side and the positive 
sign on the FM side. The resulting structure on the polar 
side is that of Eq. (|22]) , in which the nematic axis d forms 
a radial hedgehog d = r [50, 56]. This represents the 
polar point defect on the interface. On the FM side the 
spinor is similar to the singular spin vortex of Eq. ([TT]) 
with vortex lines of opposite winding in (±. Hence we 
have constructed on the FM side an approximation to a 
spin vortex that terminates to the polar monopole at the 
interface, as illustrated in Fig. SJ^d). 

Next we show that vortices can also be made to ter- 
minate at the interface. In Eq. and in Fig. Hl^b), 
a singular, singly quantized polar vortex perforates the 
interface to a coreless FM vortex when a 7r-phase kink is 
inserted in The resulting defect can be cut in half 
while still preserving the coherent interface with a con- 
tinuous order-parameter field by inserting an additional 
phase kink in ("q. This allows the vortices on different 
sides of the interface to move apart: 

/ -sin/3 \ 
C'^^ = V2e'^ COS f3 , for z > (27a) 

V2 \^ e^^^sin/^ / 

_ / sin/3 \ 
C'^* = V2e'^ COS f3 , for z < 0, (27b) 

V2 y e^^^sin/^ J 

where we may choose /3 = 37r/4 as in Eq. (|24|) . One 
possible configuration is illustrated in Fig. [5fa), where 
the singular polar vortex and a doubly quantized FM 
vortex are spatially separated and both terminate on the 
interface. 

Since the vortex lines in the individual spinor compo- 
nents terminate on the soliton planes, it is also possible 
to consider a state where a vortex exists only on one side 
of the interface, for instance, 

^ / -sin/3 \ 
= —[ V2e^^cos/3 , for z > (28a) 
V2 y e^'^sm/3 J 

[ sin/3 \ 

= 1= V^cos/3 , for z < 0. (28b) 

v2 y sin/3 I 

In addition to singly quantized vortices also half- 
quantum vortices are possible in the polar phase 
[Eq. (f2T]) ]. A singly quantized vortex can split into two 
half-quantum vortices while preserving the topology, and 
such splitting can be energetically favorable [63] . The de- 
fect configuration of perforating singly quantized vortices 
in Fig. m^a) can therefore also deform to a state in which 
a singly quantized FM vortex continuously connects to a 
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FIG. 5. (Color online) Schematic illustrations of more com- 
plicated vortex connections that can be constructed from the 
basic defect solutions illustrated in Fig. (4] (a) A singly quan- 
tized vortex in the polar phase connecting to a doubly quan- 
tized vortex on the FM side may be cut in half at the interface 
and the resulting vortices in the two regions may be moved 
apart if an additional dark soliton plane is introduced in Co- 
(b) A singly quantized polar vortex can split into two half- 
quantum vortices. Such a splitting may for example happen 
as a result of energy relaxation. This can yield a state where a 
FM vortex connects to a pair of half-quantum vortices, (c) A 
Dirac dipole can be constructed by joining the Dirac strings 
of a Dirac monopole and an antimonopole 57]. Placed on 
the interface, the dipole connects to two singly quantized vor- 
tices on the polar side, (d) Rotation may cause interface- 
crossing vortex structures to nucleate. Here four nucleated 
vortex complexes are drawn together with a singular vortex 
already present. 



pair of half-quantum vortices as illustrated in Fig. [5fb). 
We will demonstrate in the next section that this state 
does indeed appear as a consequence of energy minimiza- 
tion. 

The vortices and monopoles in the spin-1 BEG are 
made up of vortex lines and soliton planes in the in- 
dividual spinor components. The construction of the 
interface-crossing defect solutions was achieved by identi- 
fying defects in the two phases that have a similar combi- 
nation of vortex lines, and using soliton planes to achieve 
the switch between polar and FM manifolds. Dark soli- 
tons have been phase-imprinted experimentally [89. ^QOj. 
Phase-imprinting of vortex lines in a BEG by transfer- 
ring angular momentum from an electromagnetic field 
has been proposed theoretically [9?-^ and several of the 
techniques have also been realized in experiments 
99]. In spinor BEGs, coreless vortices have been pre- 
pared by adiabatic ramping of a magnetic field along the 
trap axis [7Q|-[72|. More complicated vortices and tex- 
tures could be imprinted using proposed techniques for 
creating vortex rings 
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Other interface defect configurations can be con- 
structed by combining more elementary defect connec- 
tions. For instance, in the FM phase the Dirac monopole 
can be turned inside out to form an antimonopole, in 
which the spins point radiahy inward. By joining the 
Dirac strings of a Dirac monopole and an antimonopole, 
a Dirac dipole can be constructed In Fig. Efc), 

we illustrate a state where such a dipole is placed on 
the polar-FM interface. The Dirac string forms a dou- 
bly quantized vortex line in the FM phase, connecting 
the two monopoles. Consequently, the Dirac dipole can 
form the termination points of two oppositely winding 
singly quantized vortices [see the construction that lead 
to Eq. (|25j)]. 



IV. CORE STRUCTURE OF 
INTERFACE-CROSSING DEFECTS 

A. Core deformation of interface-crossing defect 
solutions 

In the preceding section we constructed the prototype 
spinor wave functions for the interface-crossing defect so- 
lutions connecting defects in the FM and polar phases. 
We showed that such solutions can be formed by com- 
binations of elementary vortex lines and dark soliton 
(phase kink) planes that could be experimentally pre- 
pared by phase-imprinting. Here we use the constructed 
spinor wave functions for the defect configurations as ini- 
tial states for numerical studies of the defect stability. 
Provided that prototype spinors sufficiently closely ap- 
proximate the local energetic minimum configuration, the 
initial states quickly relax to the targeted defect struc- 
ture. 

By numerical simulations we can determine the en- 
ergetically preferred core structures and the energetic 
stability of the defects. In order to do so we minimize 
the free energy in the rotating frame [Eq. (j4])] by prop- 
agating the coupled Gross-Pitaevskii equations, derived 
from Eq. (0, in imaginary time using a split-step al- 
gorithm jlQl| . We assume the slightly elongated trap, 
defined by Eq. (j3j). The initial state prototype spinor 
wave functions are given in Sec. IIIIBI We choose the 
spin-independent nonlinearity cq = 2.0 x lO^hwl^, where 
/ = {h/muj)^^'^ is the transverse oscillator length. For 
^^Rb in a trap with = 27r x 10 Hz these parameters 
correspond to 10^ atoms. The spin-dependent nonlinear- 
ity, C2, is allowed to vary. 

In Ref. [63] it was demonstrated that the core of a 
singly quantized vortex in the FM phase of spin-1 EEC 
deforms by locally rotating the spin vector so that the 
vortex lines in the individual spinor components in the 
appropriate basis representation move apart. The singu- 
lar vortex line then no longer represents a vanishing atom 
density, but is occupied by atoms with zero spin magni- 
tude as in the polar phase of the spin-1 BEG. The FM 
vortex line singularity filled by atoms in the polar phase 



becomes energetically favorable by allowing a larger core 
size and a correspondingly smaller bending energy. The 
core deformation can be understood from the energet- 
ics associated with the hierarchy of characteristic length 
scales determined by the interaction strengths: the size 
of the filled core is determined by the spin healing length 
^F, Eq. ([7]), which is usually larger than the density heal- 
ing length Eq. (|6|), that sets the size of a core with 
vanishing density. 

Similarly, a singly quantized vortex in the polar phase 
was shown to lower its energy by spontaneously break- 
ing axial symmetry, splitting into a pair of singular half- 
quantum vortices. This again avoids depleting the atom 
density in the vortex core: at the location of the singu- 
larities 1(F) I = 1, and the two vortices form an extended 
core region where the order parameter is excited out of 
the polar ground-state manifold. The size of the core 
region is then enlarged to be on the order of ^f, with 
a corresponding decrease in bending energy. The over- 
all topology is preserved away from the two singularities. 
Inside the extended core region, the splitting of the singly 
quantized vortex locally deforms the nematic d field that 
describes the order parameter [see Eq. (pTj) ]. and a discli- 
nation plane, where the identification d ^ — d is made, 
appears between the vortex fines. Thus on a loop en- 
circling only one line singularity, both condensate phase 
and nematic axis wind by tt. 

It was demonstrated in Eq. (|23|) how singly quantized 
vortices in the FM and polar phases can be connected 
across the polar-FM interface, despite the fact that these 
are two topologically different defects. As the energy 
is minimized, this interface-crossing defect deforms by 
a mechanism analogous to that described above for the 
singly quantized vortices in the purely FM and polar 
BEGs. The resulting structure is shown in Fig. [6l On 
the polar side of the interface, the splitting of the singly 
quantized vortex into two half-quantum vortices is rec- 
ognized from the deformation of the nematic field, which 
shows the characteristic tt winding around each singu- 
larity, and the formation of the disclination plane. The 
order-parameter is excited out of the |(F)| = ground- 
state manifold, to reach |(F)| = 1 at the singular lines. 
The pair of half-quantum vortices connect across the in- 
terface to the singly quantized FM vortex. This, in turn, 
exhibits the local rotation of the spin vector, allowing the 
core region to fill by mixing FM and polar phases, with 
I (F) I = on the singularity. 

The relaxed interface-crossing vortex structure is thus 
recognized as that illustrated schematically in Fig. [5fb), 
and the deformation is understood in terms of the charac- 
teristic length scales set by the atom-atom interactions. 
In the purely polar or FM condensate, the core-deformed, 
singly quantized vortices are energetically stable [63] . For 
the parameter values investigated, the configuration in 
Fig. [6] ultimately decays for very long relaxation times in 
our simulations. 

We do, however, find an energetically stable deforma- 
tion of a singly quantized polar vortex connecting across 





FIG. 6. (Color online) Deformation of connecting singly quan- 
tized vortices. The green (gray) |(F)| isosurface shows the 
interface and the cores of singular defects. The polar phase 
is above the interface. Energy relaxation causes the singly 
quantized vortex on the polar side to split into a pair of half- 
quantum vortices with | (F) | > cores on the polar side, rec- 
ognizable by the characteristic tt winding of d (red/dark gray 
cones). These connect across the interface to a singly quan- 
tized FM vortex [cf. Fig. [Hb)]. A local rotation of the spin 
vector (light gray cones) allows the core region to fill with 
atoms with |(F)| < 1. (|c2| = 1.0 x lO'^hujl^ and Q = 0.20cj.) 



the interface if instead of starting from (^^^ of Eq (|23|) . 
we minimize the energy of ("^^^^ from Eq. (|24|) . This 
spinor describes a singly quantized polar vortex con- 
necting across the interface to a coreless vortex. [A 
topologically equivalent configuration can be constructed 
by allowing the polar vortex to terminate on a Dirac 
monopole, Eq. (|25]) .] Minimizing the energy leads to the 
deformation shown in Fig. [71 On the FM side of the inter- 
face, the spin structure acquires the fountain-like struc- 
ture characteristic of the coreless vortex, as shown by the 
white arrows in Fig.[7fa). Here the frequency of rotation 
determines the direction of the spin vector at the edge 
of the cloud as the angular momentum in the FM phase 
adapts to the imposed rotation. 

On the polar side of the interface we recognize the 
splitting of the singly quantized vortex into a pair of 
half-quantum vortices, identified by the tt winding of the 
nematic axis d around each vortex [Fig. El^b)], preserv- 
ing the overall topology. As before, in the core region 
the order parameter is excited out of the ground-state 
manifold, with |(F)| = 1 on the singular lines. Fig. [Tfa) 
shows how the spin texture connects smoothly across the 
interface. 

The continuity of the relaxed spinor wave function is 
further demonstrated in Fig. [51 giving a detailed picture 
of the interface region. In the relaxed state, the total 
atom density remains nonvanishing at the interface and 
varies smoothly across it. The populations of the indi- 



FIG. 7. (Color online) Minimizing the energy of Eq. ([24)) . 
corresponding to Fig. [HJb) results in a splitting of the singly 
quantized polar vortex, while the characteristic fountain-like 
spin structure of the coreless vortex is established in the FM 
part of the cloud, (a) The magnitude of the spin (|(F)| = 1 
is dark red), shows the interface after relaxation of the en- 
ergy, and the filled cores, with |(F)| = 1 at the singularity, 
of two half-quantum vortices in the polar part. White ar- 
rows show the spin vector and indicate the coreless vortex 
in the FM part, and how the spin structure connects to the 
vortex cores across the interface. This result was obtained 
using \c2\ = 2.5 x lO'^hwl^ and Q = 0.12cj. (b) The half- 
quantum vortices may be identified by the winding of the 
nematic axis d (unoriented but shown as cones to emphasize 
winding). This displays the characteristic d ^ — d winding as 
any single vortex core is encircled. The two cores are joined by 
a disclination plane. Note that away from the core region the 
original topology of the singly quantized vortex is preserved. 
Here a stronger spin-dependent nonlinearity has been used 
to get more sharply defined FM cores: |c2 I = 1.0 X lO^/kj/^ 
Q - 0.19CJ. 



vidual spinor components are also continuous across the 
interface. The continuity of the defect-carrying spinor 
wave function as it crosses between the different bro- 
ken symmetries means that the defects connect continu- 
ously across the interface. In Fig [8l the position where 
C2 changes sign is indicated by a dashed line, and |(F)| 
shown in panel (e) shows the finite width of the interface 
region after energy relaxation. 

The vortex core structures are particularly intriguing 
and complex when a vortex terminates to a point de- 
fect on the interface. We study a singular FM vor- 
tex terminating to a radial hedgehog at the interface 
(d forms a hemispherical hedgehog on the polar side), 
as depicted schematically in Fig. [H^d). Unlike the FM 
Dirac monopole, the polar monopole cannot unwind into 
simpler vortex configuration. Therefore, while the con- 
nection depicted in Fig. [Hb) is equivalent to Fig. [H^c) 
and leads to the same energy-minimizing defect configu- 
ration. Fig. [11(d) is topologically distinct from all other 
vortex connections. 

A constructed prototype spinor wave function repre- 
senting a singular FM vortex terminating to a hedgehog 
point defect is given by Eq. ([26]) . At the point defect 
singularity the atom density is zero. The density de- 
pletion at the defect core is energetically costly, and if 
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FIG. 8. (Color online) The condensate wave function varies 
smoothly across the polar-FM interface, showing that the de- 
fect states cross the interface smoothly, connecting the two 
topologies, (a)-(c) Densities in the individual spinor compo- 
nents C+5C0 and C-5 respectively. These vary smoothly and 
do not vanish simultaneously on the interface. The position 
of the interface, defined as the plane where C2 changes sign, 
is indicated by the dashed line, (d) The smooth variation 
of the nonzero atomic density across the interface shows the 
continuity of the condensate wave function, (e) (F) shown as 
color map (dark red at |(F)| = 1) and vector field (white ar- 
rows) . Note how the spin structure connects smoothly across 
the interface, in particular connecting the FM cores of the 
half-quantum vortices continuously to the fountain structure 
of the spin in the coreless vortex. The magnitude of the spin 
also shows how the interface has acquired a finite width. 



^F, Eq. ([7j), is sufficiently large in comparison with 
Eq. (j6j), the energy cost can be reduced by deforming 
the point defect into a semi-circular line defect whose 
ends attach to the interface. The resulting arch-like de- 
fect is shown in Fig. [9] together with the spin structure 
in the FM core (a) and the nematic axis away from the 
defect on the polar side of the interface (b). The defor- 
mation of the defect, schematically illustrated in Fig. IT0\ 
is local and the topological charge of the monopole is re- 
tained: away from the defect, the radial hedgehog struc- 
ture of d is preserved. ^This implies that on any closed 
loop through the arch, d must turn by tt. Consequently, 
single- valuedness of the spinor wave function requires the 
condensate phase to also turn by tt, and we infer that the 
arch-shaped line defect is a half-quantum vortex. This 
Alice arch resembles the upper hemispheric part of the 



Alice ring — a closed half-quantum vortex ring that ex- 
hibits the topological charge of a point defect over any 
surface enclosing the def ect jSq . Alice rings also appear 
in high energy physics [102] with a topological charge 
similar to the magnetic "Cheshire" charge [103]. 

The deformation of the point defect at the interface 
to an Alice arch results from a complex interplay be- 
tween the energetic considerations, topology, and the 
length scale hierarchy in the system. The two character- 
istic healing lengths determine whether the energy cost 
of forming the extended vortex core region, where the 
singularity is occupied by the atoms in the FM phase, 
offsets the energy cost of the density depletion at a point 
defect (with a vanishing density at the singularity). We 
find an abrupt transition point to an Alice arch when 
C2 < 0.5co. For larger values of 02 the point defect at 
the interface is preferred to the arch defect. The sharp 
threshold for the deformation can be understood by the 
topology. For the point defect there cannot be a sin- 
gular point with |(F)| = 1 and a spherically symmetric 
core region with |(F)| > 0. The nematic axis forms a ra- 
dial hedgehog structure and this would then imply that 
the spin vector (always orthogonal to d) would have to 
form a continuous tangent vector field for the spherically 
symmetric object which is not possible ("hairy ball" the- 
orem). Therefore the point defect has to deform to a ring 
or an arch-shaped defect at the interface before the tran- 
sition from the zero density singularity to the |(F)| = 1 
singularity with nonvanishing density is possible. The 
sharp deformation threshold also appears in the stability 
analysis of a singular point defect and the Alice ring [56| 

We note that the arch defect is unstable towards drift- 
ing out of the cloud due to the density gradient resulting 
from the harmonic trapping potential. It could poten- 
tially be stabilized using a weak pinning laser to create 
a small density depletion at the center of the trap. 

The structure of the arch defect emerging from the 
point defect may be understood by studying the individ- 
ual spinor wave function components. In the point defect 
with vanishing density at the singularity, the overlapping 
vortex lines in (± intersect with the soliton plane in ("q. 
In the prototype spinor wave function (^sv^pm addi- 
tional soliton plane is present in ("+ to account for the 
switch from the polar to the FM side of the interface. 
In Ref. [56] the deformation of the spherically symmetric 
point defect was explained by a local separation of the 
vortex lines in (± such that they no longer overlap at the 
soliton plane. Here the additional soliton plane in ("+ cuts 
the vortex line in at the interface, separating the two 
parts. The spinor components (± are shown in Fig. [91(c). 
The positions of the Alice arch and the FM vortex are 
indicated, and we find how in the polar part of the cloud, 
the separated vortex lines make up the semi-circular half- 
quantum vortex. On the FM side, the unbroken vortex 
line in and the z < part of the vortex line in form 
the spinor wave function of the FM vortex whose core is 
filled by the vortex- free Co component. 

At high rotation frequencies of the trap, we find nu- 





FIG. 9. (Color online) In Fig. [UJd) a polar point defect forms the termination of a FM spin vortex. As the energy relaxes, the 
point defect deforms into a an arch-like, singular defect line that connects to the interface at both ends, (a) An isosurface of the 
spin magnitude is shown in green (light gray). The spin magnitude rises to 1 (dark red/black) on the FM side of the interface 
(z < 0) and inside the line defect. Gray cones indicate the spin vector, (b) Away from the defect, the d vector field (red cones) 
in the polar phase retains the hedgehog structure of the original monopole. (c) In the trial wave function, Eq. ([26]), both the 
FM spin vortex and the polar monopole are formed by exactly overlapping vortex lines of opposite winding in ^^^p™ and 
^sv^pm^ Constant-density surfaces for n (red/medium gray) and n (blue/dark gray) show how the arch-shaped line 
defect is formed by deformation of these vortex lines, such that they no longer overlap close to the interface. The half-quantum 
vortex line (above the interface) and the spin vortex (below the interface) are indicated by silver and gold (light gray) spin 
isosurfaces at |(F)| = 0.9 and |(F)| = 0.5, respectively. (|c2| = 5.0 x lO^/kj/^ and Q = 0.) 




FIG. 10. (Color online). Deformation of the hedgehog point 
defect (left) on the interface into an arch-shaped line defect 
(right). Red cones represent the (unoriented) d field; the de- 
fect cores are shown as a green sphere and torus, respectively. 
The deformation creates a hole in the point-defect core, allow- 
ing it to expand into the line defect by locally deforming the 
d field. Away from the line defect, the asymptotic hedgehog 
structure is preserved. (Dashed lines indicate unperturbed 
d field for reference.) On any closed loop through the arch, 
d therefore winds to — d, identifying the defect as a half- 
quantum vortex line. The line-defect core is filled with atoms 
with 1(F) I > 0; spin vector indicated by silver arrows. 



cleation of interface-crossing defects in the energy min- 
imization. Provided that the appropriate instability for 
nucleating vortices is triggered, the emergence of defect 
configurations where a half-quantum vortex connects to 
a coreless vortex spontaneously emerges, as they lower 
the energy of the system in a sufficiently rapidly rotat- 
ing trap. An example is shown in Fig. [TTJ in which four 
interface-crossing vortices nucleate. On the polar side 
of the interface, four singular lines appear, on which 
1(F) I = 1. These vortices may be identified as half- 
quantum vortices through the winding of the nematic 




FIG. 11. (Color online) Left: Isourface of |(F)| showing 
the cores of singular defects, together with d (red/dark gray 
cones) on the polar side of the interface, and (F) (light gray 
cones) in the FM part. Four nucleated vortices are identified 
as polar half-quantum vortices connecting across the interface 
to nonsingular coreless vortices. Right: Top-down view of (F) 
in the FM part, showing the four coreless vortices. 



axis. On the FM side, the order parameter remains non- 
singular. However, the spin texture reveals that each 
of the four half-quantum vortices connects to a coreless 
vortex. 



V. CONCLUSIONS 

Interfaces between topologically distinct ground-state 
manifolds play an important role in several areas of 
physics, many of which are difficult or impossible to 
access experimentally. Here we have given a detailed 
analysis of our proposal for how topological interface 
physics of defects and textures can be studied in ultra- 
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cold atomic gases [23]. We considered the specific exam- 
ple of a spin-1 BEC with spatially separated polar and 
FM regions within a continuous condensate wave func- 
tion. For a number of representative, interface-crossing 
defect states, we provided detailed constructions of proto- 
type spinor wave functions by considering how basic vor- 
tex and monopole solutions with similar spinor structure 
represent different topological objects in the two phases. 
The resulting wave functions are built from vortex lines 
and soliton planes in the individual spinor components 
that could be phase- iniprinte d usi n^ existing experimen- 
tal techniques [45|, [H, H lOGj-ITooj . 

The energetic stability and energy-minimizing core 
structures of the interface-crossing defect configurations 
were analyzed by numerically minimizing the energy of 
the prototype wave functions. The resulting spinor states 
demonstrate how the condensate wave function smoothly 
interpolates between the two ground-state manifolds by 
locally restoring its full symmetry, thereby establishing 
a coherent interface through which defects may connect 
continuously. In particular, we demonstrated the ener- 
getically stable connection of a coreless vortex to a pair 
of half-quantum vortices, and the formation of an Alice 
arch: the deformation of a point defect at the interface 
into an arch-shaped half-quantum vortex line that pre- 
serves the topological charge. 

In order to demonstrate the basic principle of the topo- 
logical interface physics in ultracold atoms, we have con- 
centrated in this work on a relatively simple and accessi- 
ble example of defect perforation across constructed in- 
terfaces in spin-1 BECs. The interface analysis, how- 
ever, can also be applied to more complex systems, such 
as spin-2 [13, [28, 64] and spin-3 ^ loj BECs, where, 
for example, non-Abelian defects are predicted [65|, [66|. 
Other particularly promising platforms for topological in- 
terface s tudies ar e strongly correlated atoms in optical 
lattices [1 04-" 106] exhibiting also quantum phase transi- 
tions and potential analogues of exotic superconductiv- 
ity [lo| in crystal lattices. 

Moreover, the interface scheme may be used to inves- 



tigate nonequilibrium dynamical scenarios for produc- 
tion of topological defects and textures in phase tran- 
sitions. An intriguing possibility is production of topo- 
logical defects in experiments inspired by brane-inflation 
models, where brane annihilation leads to formation of 
defects [1, Q. In a FM BEC, a region of polar phase can 
be created by locally shifting the spin-dependent interac- 
tion strength. The resulting phase boundaries then form 
two-dimensional analogues of D-branes. When the inter- 
action shift is removed, the polar region collapses, simu- 
lating brane-antibrane annihilation, and resulting defects 
can be observed in the cloud. This is similar to a recent 
experiment in ^He [14], where, however the defects are 
more difficult to observe directly. 

Defect production may also result from dynamic in- 
stabilities, for example arising from superfluid counter- 
flow between the FM and polar phases. The boundary 
between two fluids moving with respect to each other 
becomes unstable if the relative velocity exceeds some 
critical value, leading to excitations on the interface. 
This phenomenon is well understood in classical fluid 
mechanics and is known as the Kelvin-Helmholtz insta- 
bility p^OTj . An analogous superfluid Kelvin-Helmholtz 
instability has been shown to occur at the interface be- 
tween superfluid ^He A and 5, providing another active 
area of research related to interfaces between different 
ordered phases . Vortices in the A phase cause coun- 
terfiow against the initially vortex- free B phase. As the 
relative velocity exceeds a critical value, vortices nucle- 
ate from the interface into the B phase. A superfluid 
Kelvin-Helmholtz instability has also been pred icted in 
phase-separated two-component BECs [ lOSl llOOj . 
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